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Abstract 

We study a model-independent parameterization of tiie vector pion form factor 
that arises from the constraints of analyticity and unitarity. Our description should 
be suitable up to ^/s ~ 1.2 GeV and allows a model-independent determination of the 
mass of the p{770) resonance, Mp = (775.1 it 0.5) MeV. We analyse the experimental 
data on r~ — > tt~tt^i't-, in this framework, and its consequences on the low-energy 
observables worked out by chiral perturbation theory. An evaluation of the two pion 
contribution to the anomalous magnetic moment of the muon, a^, and to the fine 
structure constant, a(M^), is also performed. 



PACS numbers: 13.35.Dx,13.40.Gp,13.65.+i,12.39.Fe 



1 Introduction 



The hadronic matrix elements of Quantum Chromodynamics (QCD) currents play a 
basic role in the understanding of electroweak processes at the low-energy regime (typically 
E 1 GeV). However our poor knowledge of the QCD dynamics at these energies introduces 
annoying and serious incertitudes in the description and prediction of the processes involved. 

To bypass this problem several procedures have been addressed in the literature on this 
topic. On one side there is a widespread set of models that pretend to describe, in a simplified 
way, the involved dynamics |^. While of importance to get a feeling of the entangled 
physics, the included simplifying assumptions are usually poorly justified and, sometimes, 
even inconsistent with QCD. Ad hoc parameterizations of the matrix elements have also 
been extensively used 0, ^]. The problem with this technique is that, while the description 
of data can be properly accounted for, it is not easy to work out the physics hidden in the 
parameters. 

A more promising and model-independent procedure is the use of effective actions from 
QCD. At very low energies {E <^ Mp, with Mp the mass of the p(770) resonance) the most 
important QCD feature is its chiral symmetry that is realized in chiral perturbation theory 
(xPT) [|], a perturbative quantum field theory that provides the effective action of QCD in 
terms of the lightest pseudoscalar mesons. xPT has a long and successful set of predictions 
both in strong and electroweak processes |0. At higher energies {E ~ Mp), resonance chiral 
theory is the analogous framework where the lightest resonance fields are kept as explicit 
degrees of freedom. With the addition of dynamical constraints coming from short-distance 
QCD, resonance chiral theory becomes a predictive model-independent approach to work 
with. 

One of the simplest hadronic matrix elements of a QCD current is the vector pion form 
factor Fv{s) defined through, 

{nHp)n~{p')\V;\0) = {p-p')pFy{s) , (1) 

where s = q'^ = {p + p')"^ and is the third component of the vector current associated 
with the approximate SU{3)v flavour symmetry of the QCD lagrangian. The vector pion 
form factor drives the hadronic part of both e~^e~ tt+tt^ and Ti^Ti^Ur processes in 

the isospin limit 0. There is an extensive bibliography on the study of this form factor that 
we do not review in detail here. 

At very low energies, Fv{s) has been calculated in xPT up to 0{p^) [|^, |[. A successful 
study at the p(770) energy scale has been carried out in the framework of the resonance 
chiral theory (the effective action of QCD at the resonance region) in Ref. |p. In this 
last reference the unitarity and analyticity properties of the vector pion form factor were 
implemented in order to match the low-energy result at 0{p'^) in xPT with the correct 
behaviour at the p(770) peak. The result is in excellent agreement with the data coming 

^If isospin symmetry is broken, there is a mixing between the third and eighth components of the vector 
current. The spectral functions are then sUghtly different in e~^e~ annihilation and tau decays. 



1 



from e+e" — > tt+tt" and e~7r='= — *• e~7r^ processes. This solution, that includes the p(770) 
contribution only, leaves just one free parameter, Mp, and provides a suitable description of 
Fv{s) up to y/s ~ 1 GeV. If we want to be able to extend its validity at higher energies we 
should take into account other contributions. To achieve this feature, the analyticity and 
unitarity properties of Fv{s), together with the resonance chiral theory, continue to provide 
a model-independent solution for the vector pion form factor that we analyse, in detail, in 
this article. The new solution includes two, a priori, unknown parameters in addition to Mp. 
These parameters happen to be related to the chiral low-energy observables in Ref. 0, ^j, 
the squared charged pion radius, and the O(s^) term in the chiral expansion Cy. 

In the next section we construct the vector pion form factor on grounds of its analyticity 
and unitarity relations. In Section 3 we study the experimental data on r~ 7r~7r*^z/^ with 
our solution for the pion form factor. By a fitting procedure we determine the values of Mp 
and low-energy parameters that tau decay data demand. Section 4 is devoted to analyse the 
results we have got from the fitting procedure and the consequences on the chiral observables 
of xPT. A corresponding evaluation of the two-pion contribution to the hadronic part of 
the anomalous magnetic moment of the muon and the fine structure constant is collected in 
Section 5. We present our conclusions in Section 6. 



2 Analyticity and unitarity in Fv{s) 

The vector pion form factor Fv{s) is an analytic function in the whole complex s-plane, 
but for the cut along the positive real axis, starting at the lowest threshold s = 4m^, where 
its imaginary part develops a discontinuity. This is given by the unitarity condition 

ImFvis) = \Y. I rfp„(vr+7r-|Tt|n)(n|Vj|0) , (2) 

n •' 

where \n) represents on-shell intermediate states and T"!" is the scattering operator connecting 
the intermediate state \n) to the final two-pion state. The first allowed intermediate states 
are 27r, At: and KK. To every intermediate state corresponds a branch point at the value 
of s equal to the squared sum of masses of the corresponding particles, i.e. s = (2m^)^, 
s = (4m,r)^, and so on. In the elastic region, s < 16m^, the only intermediate state 
considered in Eq. (|^) is the one with 2it, and Watson final-state theorem |10] relates the 
imaginary part of Fv{s) to the partial wave amplitude tl{s) for vrvr elastic scattering with 
angular momentum and isospin equal to one. Thus, from Eq. (|^), 



ImFvis + ie) = t}(s) = e*^i sm 6] Fyis)* , (3) 



where cr^ = y 1 — 4m^/s. As lmFv{s + ie) is a real quantity, the phase of -FV(s) must be 
(5j(s), that is, the phase-shift of the t\{s) partial wave amplitude. Therefore, 

Im Fy (s + ie) = tan (5j Re (s) . (4) 
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The analyticity and unitarity properties of Fv{s) are accomplished by demanding that the 
form factor should satisfy a n-subtracted dispersion relation in the form 



d'' ^ . ,, s" dz tan 6l{z)ReFviz) 



s 



Ms) = }2 T7TtMs)\s=0 + - / - , (5) 

k\ ds" vr Jimi z"^ z - s-ie 



where we have used Eq. (H). This integral equation has the known Omnes solution & 11 



^ / X ^ / X { s"" dz 5}{z) \ 

F.W = q„Wexp|-£__^_-,-^| . (6) 

with 

Qn{s) = exp { E ^ ^ \^Fy{s)l=o } • (7) 

Strictly speaking the solution (||) for Fv{s) is valid only below the inelastic threshold (s < 
16m^). This is because we have only included the two-pion threshold in the unitarity relation 
d^). However, the contributions from higher multiplicity intermediate states are suppressed 
by phase space and ordinary chiral counting. 

As in any subtracted dispersion relation like the one given by Eq. (^) there is an inter- 
play between the subtraction constants (polynomial part) and the dispersive integral. By 
increasing the number of subtractions (correspondingly increasing the power of z in the de- 
nominator) we pull in the low-energy part of \mFy{s) in the integrand . Then the values 
of \mFv{s) in the upper part of the integration are less important. At the same time the 
information of this high energy region shifts to the increasing number of subtraction con- 
stants that are related with the low-energy expansion of the form factor. This situation 
is reflected in the solution of the integral equation (^). If we know the 5\{s) phase-shift 
only at very low energies, an accurate evaluation of the integral in Eq. (^ would require a 
high number of subtractions. This exchange of information between high and low energies 
is, by no means, paradoxical. It is a strict consequence of the fact that, being an analytic 
function in the complex s-plane, the behaviour of Fv{s) at different energy scales is related. 
Dispersion relations embody rigorously this property. 

The 6\{s) phase-shift is rather well known, experimentally, up to ~ 2 GeV. Resonance 
chiral theory provides a model-independent analytic expression that describes properly the 
p(770) contribution H] to it : 



5\{s) = arctan 



Ml-s 



with T p{s) the hadronic off-shell p width |jT2| (see Eq. ( |A.2D in the Appendix). This result. 



that provides our definition of Mp, follows from Eq. (^) and the expression for Fv{s) obtained 
in Ref. that we collect in Appendix A. The description of data given by 5\{s) in Eq. (H) 
is accurate enough up to ii^ ~ 1 GeV for values of Mp in the ballpark of the average value 
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collected in the Review of Particle Properties [|13|. At higher energies heavier resonances 
with the same quantum numbers pop up and to get a correct description we should use the 
available experimental data from Ochs fl^ . 

We will take the result for Fv{s) in Eq. (^) with 3 subtractions. There are several reasons 
to take this case. On one side the number of subtractions is high enough to weight the low- 
energy behaviour of Sl{s) that is much well known than its high energy part. On the other 
side the number of subtraction constants, three a priori unknown parameters, is low enough 
to allow a reasonable parameterization. In fact one of the subtraction constants is provided 
by the normalization condition on the form factor, i.e. Fv{0) = 1, and there remain two 
parameters that can be related to the low-energy expansion of the form factor, (r^)y and 
Cy, as we will shortly see. 

Therefore we take as the vector pion form factor provided by analyticity and unitarity 
the expression 

Fv{s) = exp{ ais + -a2S^ + — / } . 9 

[ 2 TT J Ami Z - S - IS ) 

Since Eq. is only valid in the elastic region, we have introduced an upper cut in the 
integration, A. This cut-off has to be taken high enough not to spoil the, a priori, infinite 
interval of integration, but low enough that the integrand is well known in the interval. 
As commented above we know best 6l{s) up to -E < 2 GeV. We will take A = 2.0 GeV 
though, with three subtractions, there is a negligible difference (within the errors) between 
A = 1.5 GeV and the previous value. 

The two subtraction constants ai and a2 are related with the squared charge radius of 
the pion (r^)y and the quadratic term Cy in the low-energy expansion of the pion form 
factor 

Fy{s) = 1 + ^ {rYyS + c-ys' + 0{s') , (10) 



through the relations 



{r^)v = 6ai 



cv = ^(«2 + «?) , (11) 

that follow from the expansion of the form factor in Eq. (^) and its comparison with Eq. (p^Of) . 
We will use them to predict these observables. 



3 The mass of the p(770) resonance from a fit to r decay 
data 

The fact that Fv{s) is dominated by the p(770) vector meson up to ~ 1 GeV has been 
extensively used to get the properties of this resonance. In order to proceed, a Breit-Wigner- 
like form factor is usually introduced and fitted to the data. This procedure, however, relies 
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in a modelization of the form factor that is not necessarily consistent with QCD. Here we 
propose a thorough model-independent determination of the mass of the p(770) resonance, 
Mp, defined by Eq. (|). 

Fv{s) endows the hadronic dynamics in the ti^tt^Ut decay and the e'^e' tt+tt^ 

process. The experimental data from this last source [|1^, has been available for long 
time and deeply analysed. The decay r~ — 7r~7r°z/^ has recently been measured accurately, 



in the energy region of our interest, by three experimental groups : ALEPH [|17|, CLEO-II 



m and OPAL |9|. We take Fv{s), as given by Eq. (D, to fit the ALEPH set of data. 



An appropriate study of the form factor requires a proper description of the Sl{s) phase- 
shift in the integration interval. As we are working with 3 subtractions the main contribution 
to the integration in Eq. (|^) comes from the low-energy region of the phase-shift. However if 
we wish to consider -FV(s) around y/s ~ 1 GeV the cut-off A should be not lower than, let us 
say, a/s ~ L5GeV, as we commented previously. Therefore we require a precise description 
of Sl{s) in this energy region. We achieve this through the following procedure : Sl{s) given 
by Eq. (||) provides an implementation up to ^JSmatch = Mp] hence for Mp < y/s ^1.5 GeV 
(higher values of a/s being unimportant because the three subtractions performed) we include 



the Ochs set of data fl^. As a result we come out with a description of 5{(s), in the region 
of interest, that contains all the necessary physics input. 

However there are still contributions to the form factor in Eq. (P) that are not taken into 
account with Ochs data. These are those of coupled channels that open at the KK threshold 



[20|. Therefore, in order to have a conservative determination of the observables, we choose 



to fit ALEPH data in the range 0.32 GeV < a/s < 1.1 GeV where we have a thorough control 



of the contributions. The fitting procedure is carried out with the MINUIT package pT |. 
We find 

Mp = (775.13 ± 0.02 )MeV , 
ai = (1.84 ± 0.02) GeV~^ , 
as = (4.18 ± 0.05) GeV"^ , 
xVd.o.f. = 33.8/21 . (12) 

Though the x^/d.o.f. value found can be considered reasonable it is necessary to notice 
that 80% of comes from just three points ^. Errors in Eq. (|12D , given by the MINUIT 
program, are to be taken with care. They do not include those that come from the choices 
we have made in our approach: the energy range to be fitted, number of subtractions, upper 
cut of integration A and the matching point, ^JSmatchi between Ochs data and Eq. (^. 
We estimate the final errors by exploring the stability of the results with two and four 
subtractions, varying the cut-off from A = 1.5 GeV to A = 2.0 GeV, extending the fitted 
energy range up to ^/s ~ 1.6 GeV and shifting ^Smatch within the Ochs data errors. Hence 
we conclude the figures 



^One of them at ^fs ~ 0.70 GeV and the other two around y/s ~ 0.85 GeV. 
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Mp = (775.1 ± 0.5)MeV , 
ai = (1.84 ± 0.05) GeV-2 , 

a2 = (4.2 ± 0.2)GeV-^ . (13) 

The parameters ai and 0:2 turn out to be highly anti-correlated. This procedure provides a 
mass for the p(770) resonance roughly 5 standard deviations higher than the Particle Data 
Group new average |T^ that is Mp = (769.3 ±0.8) MeV but consistent with their average 
from T decays and e+e~ processes, Mp = (776.0 ± 0.9) MeV. 




Vs (GeV) 



Figure 1: Comparison of the result of the fit to ALEPH data with the experimental ALEPH 
IT1\] and CLEO-II / |7^/ data on Fv{s) from tt^tt^z/t- in the p(770) energy region. The 



result of Ref. for Mp = 775 MeV is also shown. Up to i/i ~ 0.8 GeV both curves are 
almost indistinguishable. 

In Fig. 1 we compare the experimental data with our prescription. We also include the 
parameter-free prediction (1 subtraction only) of Ref. that, for completeness, we recall in 
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s/V(lsl) (GeV) 



Figure 2: Comparison of the result of our fit with the experimental data on Fv{s) from 
e~^e~ tt^tt" (time-like) / (J5|/ and e^n"^ e^n^ (space-like) JT^/ . The result of Ref. 
(Mp = 775 MeV) is also shown. In the region —0.4 GeV < s/ ^J\s\ ^ 0.8 GeV both curves are 
almost indistinguishable. 

Appendix A. It can be seen that our fit gives a good description of data up to ii^ ~ 1.3 GeV. 
Experimental data (in spite of tlie big errors in the higher energy region) seem to have a 
determinate structure (mild shoulder) around E ~ 1.3 GeV. This could be due to a heavier 
p-like resonance as the p(1450). Our solution takes into account this possibility though, 
because Ochs data embody these resonances up to ~ 1.5 GeV. 

We can compare the results of our fit to tau decay data with the experimental results 
coming from e~^e~ — >■ vr+vr" (time-like) [|I3| and e~7r^ e~7T^ (space-like) |TB[ processes. 
In Fig. 2 we show these sets of data together with the same curves of Fig. 1. We conclude 
that the agreement of our fit with data is good within the errors. Notice that e~^e~ — > tt+tt^ 
data has a contribution from 0^(782) that translates into a slight deformation on the right- 
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{r% (GeV-^) 


cl (GeV-^) 


Our fit 


11.04 ±0.30 


3.79 ±0.04 




11.22 ±0.41 


3.85 ±0.60 



Table 1: Low-energy observables of the vector pion form factor up to the quadratic term. 
We give the results for our fit and the 0{p^) xPT analysis of Ref.^. 

hand side of the p(770) peak. This is due to a small / = component contributing to the 
27r spectral function in e^e^ tt+tt^. This contribution does not appear in the isovector 
spectral function from -k^-k^Vt which we are describing. 

4 The low— energy observables 

At E < 2m,r the vector pion form factor satisfies a low-energy expansion given by 
Eq. ([T0|). Up to the quadratic term in s we have, therefore, two low-energy observables, the 
squared charge radius of the pion, (r'^)y, and the quadratic term Cy, that are related with 
the parameters ai and 0:2 of the form factor (^) as given in Eq. (11). 

(r^)y and Cy have recently been determined at 0{p^) in xPT While chiral symmetry 
constraints successfully provide the chiral logarithms, it remains an incertitude in the poly- 
nomial part that involves counterterms not predicted by the chiral framework. Therefore 
it is not possible to give a plain prediction for these observables. The authors of Ref. ^ 
performed, by including properly the chiral logarithms, a fit of the pion form factor, as given 
by 0{p^) xPT, to the data from t" ti~ti^Vt-, e'^e~ tt^tt~ and e~n^ — >■ e'n^ in the 
low-energy region {E < 0.5 GeV). Our procedure provides the low-energy observables from 
a fit to a larger energy interval in the time-like region. In Table 1 we compare our figures 
with those of Ref. 0. As can be seen the results compare very well but the errors to the 
observables provided by our procedure are smaller (noticeably in Cy). 

As commented above the predictability of xPT at 0{p^) is spoiled because chiral sym- 
metry does not provide information on the finite part of the counterterms in the results 
of (r^)y and Cy. Two combinations of 0{p^) counterterms, ry^(Mp) and ry^{Mp), one on 
each observable, have to be considered. In order to predict these terms one has to rely in 
modelizations or dynamical assumptions like Vector Meson Dominance (VMD). This last 
resource was employed in Ref. ^ to evaluate the vector resonance contributions ry^(Mp), 
ry^{Mp) that is the dominant piece by far. 

Numerically the 0{p^) xPT expressions relating the low-energy observables with the 
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r|;^(M,) X 10^ 


r|;^(M,) X 104 


Our fit 


-0.79 ±0.19 


1.46 ±0.03 




-0.68 ±0.26 


1.50 ±0.44 


VMD 


-0.25 


2.6 



Table 2: Combination of 0{p^) counterterms appearing in the xPT evaluation of and 
Cy . We give the predictions from our fit and the ones from the chiral fit and the VMD result 
ofRef.M- 



polynomial terms are 



12.312 ± 1603.4 r|;^(Mp) J (GeV"') 
1.787 ± 13718.7 r^^(Mp)j (GeV"^) 



(14) 



Within VMD these counterterms are obtained by integrating out vector resonances using 
the resonance chiral theory framework 0. They have been worked out, within the Proca 
formalism, in Ref. P[ with the results 



Vi 



V2 



2^2— Ay fyf^ 

■9vf\ 



V 



Ml 



(15) 



obtained by integrating the lightest octet of vector resonances of mass My. The couplings 
fv, gv and f^ can be phenomenologically obtained from p — > e+e~, p tt'^tt" and (p KK 
with the results: fy = 0.20, gy = 0.09 and f^ = —0.03, and, therefore, giving values for 
Vy, that we collect in Table 2. We compare these VMD results with the ones obtained from 
our fit and the ones provided by the 0{p^) xPT fit. We notice that the result of VMD 
seems to undervalue |ry^(Mp)| and overestimates ry^{Mp). As can be seen from Eq. ([I^) 
this difference would affect most the value of Cy. It has to be observed though that, on 
one side, to extract ry^[Mp) from (r^)y in Eq. ([I^ ) a strong cancellation driven by the term 



^For a complete discussion see Ref. B. We take for 



their Set I possibility. Our numbers differ 



slightly from the ones given in that reference because we use for the pion decay constant = 92.4 MeV 
instead of ~ 93.2 MeV. We neglect the small local contribution from pseudoscalars. 
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((r^)y — 12.312) is involved and therefore it is very sensitive to the value of the squared 
charge radius of the pion (this problem does not arise in the ry^{Mp) case); on the other 
side, VMD can only offer a rough estimate because, at this order, heavier resonances could 
also give a noticeable contribution while the VMD result only includes the lightest octet of 
vector mesons. By neglecting these heavier states we could invert the procedure and use 
our fit to predict the products of couplings f^fv and gvfv from Eq. (^). We obtain, for 
example, fx/ov = (—1-9 ± 0.6) far from the phenomenological value fx/9v — —0.33. It 
looks as if the role of heavier resonances is crucial in order to describe 0{p^) vector driven 
contributions in xPT. 



5 Two— pion contribution to the muon {g — 2) and to 

a(M| 



The hadronic contribution to the anomalous magnetic moment = {g^ — 2)/2 of the 
muon is the main source of incertitude in its theoretical prediction. Its leading part comes 
from the photon vacuum polarization insertion into the electromagnetic vertex of the muon. 
It gives |2^ : 

al^'^{vac.pol.) = (692.4 ± 6.2) X 10"^° . (16) 

This contribution can be evaluated in terms of the experimental hadronic total cross-section 
a{e~^e~ —>■ hadrons), where e"^e^ tt+tt^ is, by far, the dominant part at low energies. The 
bulk, both of the central value (~ 75%) and the error (~ 80%), of a^"''^ in Eq. (|T6|) comes 
from this vrvr intermediate state p3| . 

The relevant dispersion integral to evaluate this contribution is (up to two loops) 



'a{0)mA ds 



3 TT / J Ami S 



/ -R^^{s)K{s) 

J 4ml S 



(17) 



4 TT a ' 



where the K{s) function is given in Ref. [EB]. In terms of Fv{s) we have 



where we have introduced a cut-off A as the upper limit of integration. As K{s) grows mildly 
at high values of s, the integration in aj^^ in Eq. ([T^) is dominated by the very low-energy 
region that gives the main contribution. 

The hadronic contribution to the shift in the fine structure constant Aa{s), defined 
through a{s) = a(0)/(l — Aa{s)), can be evaluated from e~^e~ —>■ hadrons data by using 
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A (GeV) 


X 10^° 


Aa(M|)U^xl04 


1.0 


505 ±6 


33.8 ±0.4 


1.1 


511 ±6 


34.7 ±0.5 


1.2 


514 ±6 


35.1 ±0.5 


1.3 


516 ±6 


35.4 ±0.5 



Table 3: Values of al"" and Aa(M|)|^^ given by our fit to ALEPH r decay data, in the whole 
energy range (0.32 GeV < a/s < 1.6 GeV), for different values of the A cut-off. 



a dispersion relation together with the optical theorem p6[. The last estimation has been 
worked out in Ref. |^ giving 

A«(5)(M|)U,, = (276.3 ± 1.6) x 10^^ , (19) 

where the superscript indicates that only the 5 lightest quark flavours have been considered. 
The TTTT contribution can be accounted for by 



12 vr Jiral s{s-Mz) 

where, once more, we have introduced a cut-off A as the upper limit of integration in order 
to control the good description of the integrand. Contrarily to what happens in the a^^ case, 
from Eq. (|20|) we see that the integrand is not so dominated by the low-energy region and, 
therefore, higher energy contributions are relevant to evaluate Aa(M|)|7r7r. In addition, and 
as we will see, the vrvr contribution to /S.a{Mz)\had in this energy region is just a modest 10% 
of the full value (0). 

The study on the vector form factor of the pion that we have carried out allows us to put 
forward a prediction for both a"^ and Aa(M|)|7r^ that we work out as follows. The flt to 
ALEPH data that gave our results in Eq. ([T3|) was limited to ^/s < 1.1 GeV. As commented 
there we took this region because we have a thorough control of the physics involved within. 
At higher energies new physics input, unaccounted for, appears. As a result, in Fig. 1 it 
can be seen that our flt misses barely the data above i/i ~ 1.2 GeV, well outside the fltted 
region. The computation of the integrals in a'"^ (|T8|) and AQ;(M|)|,r7r (^0) requires a good 
knowledge of Fv(s) up to s ~ A^, therefore if we wish to reach A ~ 1.3 GeV we would need 
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a better description of data than the one given with the parameters in Eq. ([131) . To achieve 
this feature we fix Mp = 775.1 MeV, as concluded in Eq. (|13|), and leave ai, 02 as free 
parameters. Then we fit the ALEPH data in the whole range 0.32 GeV < ^/s < 1.6 GeV. By 
studying, as above, the stability of the fitted parameters against variations in the number 
of subtractions, the upper limit A, and the matching point y^Smatch, we conclude the values 
al = (1.83 ± 0.03) GeV-^ 02 = (4.28 ± 0.08) GeV-^ consistent with the solution of the 
restricted fit (|13D but with smaller errors. The tildes on ai and 02 are meant to prevent 
their use in Eq. ([Til). We emphasize that al and 02 are not proper physical values of the 
ai, a2 parameters because we have fitted a region of experimental data that is not properly 
implemented theoretically. However the above values of 5i and 02 describe well data up 
to -y/s ~ 1.3 GeV and, therefore, are useful to evaluate the integrals in ajl" and Aa(M|)|7r7r 
with smaller errors. The values we get are collected in Table 3. 

It has to be noticed that our errors are similar to those obtained in recent estimations 



||23|| , though the results in this reference were obtained from a combination of e~^e~ tttt 
and n^-K^Vr decay data while our results come from a fit to this last process up to 

^/s ^1.6 GeV. An improvement on our errors would require an analysis of the pion vector 
form factor with a more complete set of data, combining e+e" tttt and r~ — 7r~'7r°z/T- 
pro cesses. 



6 Conclusions 

To gain access to the resonance properties, from experimental data, a correct definition 
of those properties has to be theoretically implemented. The use of modelizations, though 
sometimes unavoidable, can spoil seriously the conclusions obtained from data. In this article 
we have studied the vector pion form factor Fv{s) within a model-independent approach. 
We have introduced a parameterization of the form factor provided by the all-important 
properties of its analyticity and unitarity relations. This last construction relates Fv{s) to 
the Sl{s) phase-shift of elastic tttt scattering. 

To proceed we have included the Sl{s) phase-shift (up to ~ 1.5 GeV) with a model- 
independent parameterization, provided by the resonance chiral theory and experimental 
data. Our form factor depends on two, a priori unknown, subtraction constants and the 
p(770) mass. We have fitted ALEPH data on tt^tt^Uj. to the form factor for E < 

1.1 GeV and we obtain Mp = (775.1 ± 0.5) MeV. Our result for Mp is bigger than the new 
average of the Review of Particle Properties [|13] but very much consistent with their average 
from r decays and e"'"e~ annihilation processes. The predictions given by our results on the 
low-energy observables worked out in xPT, (r^)y and Cy have also been computed. We find 
good agreement with the results from the fit in xPT though our errors are smaller. It is 
necessary to notice, though, that when these figures are worked on to determine local chiral 
0{p^) counterterms, the values we get are not consistent with those obtained, through VMD, 
from resonance chiral theory by integrating out the lightest octet of vector resonances. As a 
conclusion it seems that room is left for the contribution of heavier resonances. 
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Finally we have evaluated the tttt contribution to the anomalous magnetic moment of 
the muon, a^'^, and the shift of the fine structure constant Aa(M|)|7r7r. An improvement in 
the theoretical errors of these quantities would be achieved with a more complete analysis 
of the available data. 

We have shown how it is possible to extract model-independent information of resonances 
from experimental data by exploiting general properties of form factors, such as unitarity 
and analyticity. When combined with the resonance chiral theory, the effective action of 
QCD at the lightest resonance region, these properties provide a compelling framework for 
the study of form factors. 
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Appendix A 

A theoretical construction of the vector form factor of the pion was performed in Ref. 
by matching the O(p^) xPT result (valid at -C Mp) with the prescription provided by the 
resonance chiral theory. The procedure also took into account the analyticity and unitarity 
properties of Fy{s). The result only includes the contribution of the p(770) resonance and 
gives an excellent description of data up to ~ 1 GeV with just one parameter, Mp. We 
have compared this prescription with ours in Figures 1 and 2. 

For completeness we recall here the result of Ref. P] : 



Fvis) 



Mt 



—s 



Mt - 



s - iMpVp{s) 



exp 



ReA(ml/s,ml/Ml) + 



]-ReA {m\/s,m\/M^ 



(A.l) 



where ^ p{s) is the hadronic off-shell width of the p(770) resonance 

\aldis-Aml) + -aldis-Aml) 



(A.2) 



and 



A(ml/s,ml/^^^) = In (^] + 8^ - U In f ^ 



(A.3) 



with ap = Jl — Arrip / s. 
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